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TECHNICAL NOTE 2736 


TWO -DIMENSIONAL SHEAR FLOW IN A 90° EIBCW 
By James J. Kfamer and John D- Stanltz 


SUMMARY 

As part .of an approach to a better understanding of the motion of 
real fluids in flow machinery, two-dimensional, incompressible, non- 
vis cous shear flows in a 90° elbow have been investigated. Solutions ‘ 
are presented for linear and sinusoidal velocity distributions across 
the inlet of the elbow. The solutions with linear inlet velocity dis- 
tributions indicate that as the negative vorticity of the flow Increased: 
(1) the static-pressure drop through the elbow decreased, (2) the local 
deceleration along the outer channel walls increased, and (3) the mag- 
nitude of the velocities on the channel walls changed greatly, but the 
local pressure coefficient rose only gradually and the difference in 
pressure coefficient at corresponding points on the two walls was prac- 
tically unchanged. In the case of a sinusoidal inlet velocity distri- 
bution local decelerations occurred on both walls . 


’ INTRODUCTION 

As part of an approach to a better understanding of the motion of 
real fluids in flow machinery, two-dimensional, incompressible, non- 
viscous shear flows in a 90° elbow have been investigated at the NACA 
Lewis laboratory. For real, viscous fluids the velocity distribution 
upstream of the elbow is nonuniform and therefore the fluid motion is 
rotational. Such rotational, or shear, flows can develop both normal to 
and in the plane of the elbow. If the shear flow develops normal to the 
plane, so that the vorticity vectors are parallel to the plane, three- 
dimensional secondary flows develop in the elbow (references 1 to 3) • 

If, however, the shear flow develops in the plane of the elbow, so that 
the vorticity vectors are normal to the plane, the shear flow remains 
two-dimensional and in the plane of the elbow. It is this latter two- 
dimensional shear flow that is investigated in this report. The purpose 
of the investigation is to determine the effects of this type of shear 
flow on the velocity and pressure distributions in a 90° elbow. 

Two-dimensional potential solutions for irrotational flow in chan- 
nels are well known (references 4 and 5, for example) , and solutions for 
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two-dimensional shear flow about isolated todies have "been obtained 
(references 6 and 7, for example), but no solutions are known for two- 
dimensional shear flow in channels • In this report such solutions are 
presented for incompressible flow in a 90° elbow. 

Although the shear flow is considered to have been generated by 
viscous forces acting on the fluid upstream of the region of tur ning , it 
is assumed that in the immediate vicinity of the region of turning, 
where the solutions are obtained, the viscous forces are small and can 
be neglected. Solutions are obtained for two types of velocity distri- 
bution upstream of the region of turning; linear and sinusoidal distri- 
butions. streamlines and the velocity and pressure distributions of 

1 these solutions are compared with those of the potential solution from 
reference 8. 


METHOD OF ANALYSIS 


The differential equation for the distribution of shear flow in any 
two-dimensional channel is developed from Bernoulli ' s equation, and from 
the equations of motion and continuity. The flow is assumed to be 
incompressible and nonviscous. 


Equations of motion. 


- The equations of motion are 
du 


_ g dp 
p <3x 


__ Su 

u £ + v 3y 




+ V 


dv 

3y 


(la) 

(lb) 


where- g is the acceleration due to gravity; p, the fluid weight 
density; p, the static pressure; x and y, Cartesian coordinates; 
and u and v, the velocity components in the x- and y-directions, 
respectively. All symbols are defined in appendix A. All variables are 
made dimensionless by expressing them as ratios of characteristic quan- 
tities in the following manner; 


P 


P^,av } 


ns — *j _ 


- (lc) 


where q e ay . is the arit hme tic average of the resultant velocity q 
across the channel exit, downstream at infinity. 
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Qe,av 


X - — 




where w e is the width of the channel at the exit. 


(Id) 


(q~ a y) 2 

After division hy — , equations (la) and (lb) became in 


dimensionless form. 


&P 

3x 


2a 1 + 2T H 


■Sp __ dv . w &v 

" 2U 3k + 27 m 


BernO'dlli 1 s equation . - Bernoulli's equation states 

E + #.h 

P 2g 


(2a) 

(Zb) 

(3) 


where h is the Bernoulli constant (constant along a streamline), 
Bernoulli constant £s made dimensionless by expressing it as a ratio 


of h to 


_ 2g _ 


and the velocity is expressed nondimens ionally as 


a ratio of q. to q e av . Thus, 

H 


"(^e.av) 

_ 2g _ 


Q 




(3a) 


(3b) 


av 


, „ (q*. av ) 2 

Equation (3), having been divided by — , 

sional form, ° 


becomes, in nondimen- 


P + Q 2 « H 


(*) 
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Continuity equation . - The continuity equation In dimensionless 
form becomes 


dU , by 
3X + 3? 


0 


(5) 


Equation (5) Is satisfied by a stream function ¥, which Is defined by 
the following relations 


dur 

3x 


- V 



- U 


(6a) 


Also, from equation (6a) , In the direction of If 

If" Q (6b ) 

where If Is distance ' (dimensionless, expressed as ratio of w e ) along 
the outer normal to a streamline . - 


Differential equation . - From equation (4) 
U 2 + V 2 ■ 


with Q 2 


written as 


I + 20 i - ^ S 


ax 


ap «TT Su Orr by dH 

aY + 2U 5Y + 2v Sy“aY 


From equations (2a) and (7a) 



and, similarly, from equations (2b) and (7b) 



(7a) 

(7b) 

(8a) 


(6b) 


Because H is a function of ¥ only, 
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and 


dH Si? d E _ dE 

Sr “ Sr dT? “ u a? 


(9b) 


From equations (6a), (8a), and (9a) 


a 2 ^ a 2 !? i ds 


( 10 ) 


Equation (10) is the differential equation for the distribution of the 
stream function in a two-dimensional channel with shear flow. Equa- 
tion (10) can also be derived from equations (6a), (8b), and (9b). 


Vortlcity . - The varticity £ of the shear flow is defined by 


Sv Su 

£ “ ai - ar . 


(ii> 


which from equation (6a) becomes 


s 




Thus v -from equation (10) the vortlcity 1 b equal to - i 2L which is 
constant along streamlines so that 2 


, 1 dH 1 f6E\ 

5 ' 2 df " - W ) ± 

But, because the flow is parallel at the inlet. 


(Ha) 



and from equation (4) with P.j_ constant 

(§1 ■ (S) ± 

so that, after introduction of equation (6b), equation (lla) becomes 



i 


(lib) 
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Therefore, the vortieity at a point on a streamline anywhere In the 
flow field can be determined by equation (lib) from the velocity gra- 
dient normal to the streamline at the inlet. 


METHOD CEF SOLUTION 

The differential equation (10) was solved by relaxation methods for 
the examples of this report. In order to facilitate the solutions, the 
flow field in the XY -plane was transformed onto a field the coordinates 
of which are the streamlines and velocity potential lines for irrota- 
tional flow. 


Transformation of , coordinates , 
function and velocity potential, tj 
by 


- For irrotational flow the stream 
and l respectively, are defined 



(12a) 


where the subscript o Indicates irrotational flow. Also, -from equa- 
tion ( 12a) , in the direction of N 0 and S D * 



(12b) 


where Qq is the resultant velocity fear irrotational flow and where S 
is the distance (dimensionless, expressed as ratio of w e ) along stream- 
lines. From equation (12b) , If T) is zero along the right channel wall 
when faced in the direction of flow, rj is 1.0 along the left wall 
(because Qq and are 1.0 at the exit of the channel, downstream 
at Infinity). Thus, because £ extends from minus to plus infinity, 
the flow field in the £t] - plane is an infinite strip of \mit width. 

Differential equation . - The differential equation (10) transformed 
onto t!he gT] -plane' becomes (appendix B) 



Velocity components . - The velocity Q and 
and. T] Sirections, respectively, 
V by (appendix C) 


Q_ in the | i 
distribution of 


- vS? 




its components Q| and 
are obtained front the 


(14a) 



HACA TN 2736 


7 


where 


^ If 
Q n “ If 


(Mb) 

(14c) 


The velocity components U and V are related to Qg and by 
(appendix C and fig. l) 


U 


z° B e o - % — “O 


sin 8, 


V b Q| sin S Q .+ cos 0 Q 

Jf 

where 0 O is the local direction for irrotational flow (fig. 1) . 


(15) 


Bo undar y conditions . - The stream function W is constant winng 
each of the channel walls . If f is zero along the right wall (t]«* o), 
then W is 1.0 along the left wall (tjbI.O), because from equation (6b) 
and from the definition of Q 


m rn 



q air » i.o 


across the channel at the exit downstream at infinity. 


In theory the flow field extends to ±eo in the direction of £. In 
practice, however, the relaxation solutions converge to essentially uni- 
form flow conditions within a reasonably short distance from the region 
of turning. 


Relaxation solutions . - The numerical examples were solved in the 
plane ty relaxation methods (references 4 and 9) . A square grid with 
a spacing of l/Q was used. The first example was solved using a three- 
point system for computing the derivatives, and the value of UF at each 
grid point was relaxed to a unit change in the fifth place. i i t Ha second 
example was solved using a five-point system for computing the derivatives 
and the value of I at each grid point was relaxed to a unit in 

the fourth place. An investigation of the errors involved in the two 
methods of approximation to the derivatives showed comparable accuracy 
for the two solutions. 
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NUMERICAL EXAMPLES ~ 

The channel investigated by the tvo numerical examples in this 
report is the 90° elbow {89.36°) of reference 8. The elbow was designed 
in reference 8 for incompressible, potential flow with a prescribed 
velocity distribution that has no local decelerations along the channel 
walls. The transformed coordinates £,t) and the values of Qq 
required for the solution of equation (13) are tabulated in reference 8. ' 


Numerical Example. I" 

Inlet velocity distribution . - In examp lei a linear distribution 
of velocity across tin- inlet was prescribed such that the velocity Q 
varied from 0.25 on the inner wall (wall with smaller radius of curvature 
and at which t] is 0) to 0.75 on the outer wall (rj equals 1.0) • 

(See fig. 2.) Thus, 

Q ± - 0.25 + 0.50 Ti ± (16a) 

» 

For this inlet velocity distribution the. distribution of H is 
(appendix D) 


H- % +| _ (16b) 

where Hq is the arbitrary value of H along the inner wall. Therefore, 
the vorticity of the shear flow becomes 

6 " “ i S = " T ^ 16c ) 

so that for a linear inlet velocity distribution the vorticity is every- 
where constant (and negative) in the channel. Also, for the linear inlet , 
velocity distribution given by equation (16a) it can be shown that 
(appendix D) 


\ - 0.5 (tj ± + tj ± 2 ) • (l6d) 

¥ e » 0.875 T]g + 0.125 Tjg' 2 (I6e) 

.and 

Qg - 0.875 + 0.25 T] e (l6f) 

The exit velocity distribution is compared with the inlet distribution 
in figure 2. 
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Streamlines . - In figure 3 the streamlines for the shear flow of 
e xam p le I are presented in the physical XT-plane with the |,ti grid 
system, superposed. The spacing of the streamlines decreases linearly at 
hoth inlet and exit, indicating a linearly increasing variation in 
velocity at "both positions, because the spa cing of the streamlines is 
inversely proportional to the velocity. The velocity is more nearly 
u nif orm, at the exit and, therefore, at the exi t the streamlines W of 
the shear flow are closer to the streamlines tj of the potential flow 
than at the inlet. 

Velocity distribution . - Lines of constant resultant velocity Q 
are presented in figure 4.* The velocity component normal to 

potential flow streamlines is shown in figure 5* This velocity camponertt 
is relatively smal ,1 throughout the channel and disappears upstream n-n ^ 
downstream of the region of turning and along the channel walls. Two 
points of relative maximum occur. These points are explained as follows. 
From figure 1 


« Q sin (0 - 0 O ) (17) 

' f 

so that relative maxi mums for occur when the product of Q a.n ri 

sin (0 - 0 O ) Is maximum. The first martmim occurs in a region 
(tj ts 0.375, J ~ 0.625) where (0 - 0 O ) Is large (fig. 3) and the sec- 
ond maximum occurs In a region (tj ss 0.5, g S2.75) where Q Is large. 


The velocities along the inner and outer walls are plotted against 
{ in figure 6 -for example I and for the potential flow. The gradients 
of the velocity distributions for the two types of flow are similar. 
However, the magnitudes of the velocity* are different. On the outer wall 
the velocities for example I are higher than for the potential flow 
solution, whereas on the inner wall the velocities for examp le I are 
lower. These velocity distributions result from the prescribed inle t 
velocity distribution which, as indicated by equation (lib), prescribes 
the vorticity throughout the channel. Thus the circulation aro und any 
closed path, which circulation is zero for potential flow, is equal to 
£A, where A is the area of the region enclosed by the path. This 
expression for the circulation comes directly from Stokes’ theorem 
equations (5) and (11) • Although no local decelerations occur wi ring the 
channel walls for the potential flow, a slight deceleration occurs w.inn g 
the outer wall for the shear flow, example I. 


The pressure coefficient along the inner w nfl outer walls for 
example I and for the potential flow are presented in figure 7. The 
pressure coefficient is defined as 


P - P, 


i PC’le.av) | 

L 2g 


- P - P 4 


(18) 
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The pressure coefficient for example I is greater along both the inner 
and outer walls than along the corresponding "walls for potential flow, 
hut the difference in pressure coefficients on the two walls at corre- 
sponding values of | is about the same for both types of flow, indi- 
cating that the turning angle of the mean flow is about the same* It is 
interesting *to note -that for example I in the vicinity of g equals 4.5 
the pressure along the outer wall is less than along the inner wall, 
thus indicating an overturning of the average flow just ahead, of this 
region. The exit pressure coefficient is greater for example I than for 
the potential flow solution, indicating a smaller pressure drop through 
the elbow for shear flow. Physically, this smaller pressure drop for 
example I is explained by the fact that, because of the nonuniform. inlet 
velocity for example I, the average value of Hi (equation (4)) must be 
greater for example I than for the potential flow solution. But the 
linear exit velocity distribution for example I is more nearly like that 
for the potential flow solution (fig. 2); and, therefore, because the 
average value of H is higher for example I, P e is also higher. 


Numerical Example II~ 


Tnlet velocity distribution . - In example U a sinusoidal distri- _ 
bution of velocity across the inlet was prescribed (fig. 8) such that 
the velocity Q is zero on both walls and the arithmetic average value 
of Q is 0.5 (same as example I) . Thus, 

Qi » ^ sin jtt] (19a) 


This sinusoidal variation in velocity approximates the parabolic distri- 
bution for fully developed laminar flow in a straight duct with parallel 
walls. For the inlet velocity distribution given by equation (19a) the 
distribution of H is (appendix E) 

H - Ho + —Or - v 2 ) ( 10 b) 

Therefore, the vorticity of the shear flow becomes 


1 dH 

w ---- H 

2 dY 



(19c) 


so that for a sinusoidal inlet velocity distribution the vorticity is a 
linear function of W- Also, for the sinusoidal inlet velocity distri- 
bution given by equation (19a) it can be shown that (appendix E) 


¥ i “ \ ( 1 " cos ^i^ 


(I9d) 
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y e “ \ C 1 " 008 fie + 8111 fie) (19e) 

and 

% “ f (sin fl e + cos |^ e ) (I9f) 

The exit velocity distribution is compared with the inlet distribution 
in figure 8. 

Streamlines . - The streamlines of example H are presented in fig- 
ure 9. An eddy has formed on the outer wall and extends from g B -7.25 
to { ~ 2.0. (in theory this eddy extends upstream, to mi nim infinity 
but becomes extremely narrow so that within the accuracy of re laxa tion 
methods it disappears at £ ~ -7 • 25 . ) This eddy occurs because, as in 
example I, the velocity decelerates along the outer wall. Because the 
velocity was initially zero at this wall, the flow therefore reverses 
Itself and an eddy is formed. Three stagnation points occur for the 
eddy: two on the outer wall (one downstream at infinity, theoretically, 

and the other at £ approximately equal to 2.0) and a third in the 
center of the eddy nea r £ ~ -0.375. The maxi rm im velocities occur npar 
the center of the channel as Indicated by the narrow spacing of the 
streamlines. 


Velocity distribution . - Lines of constant resultant velocity Q 
are presented in figure 10. The velocity component normal to 

potential flow streamlines is shown in figure 11. r Fhl« distribution of 
Qtj appears similar to that for example I (fig. 5) . However, the mag- 
nitudes of 0^ are greater -in example U and both positive and negative 
values exist. A minimum point occurs at { " 0.375, tj Z 0.250 and a 
m a ximu m at | « 2.750, tj Z 0.625. As for example I, these points occur 
for reasons already discussed in connection with equation (17) . This 
velocity component c h a n ges sign because the angle (d - 0 O ) (see fig. 1) 
changes sign as Indicated in figure 9 where the streamlines for the shear 
flow are inclined toward the inner wall with respect to the potential 
streamlines (constant tj) in the first half, of the elbow and are in 
Inclined toward the outer wall in the second half. A. 1 irip- of zero Co- 
occurs along which the streamlines for shear flow are parallel to the 
streamlines for potential flow. 


The velocities along the inner and outer, walls are plotted against 
| in figure 12 for e x a mp le H and for the potential flow solution. Par 
the shear flow solution, the velocities on the inner and outer walls con- 
verge to zero upstream at infinity. As for example I, the flow deceler- 
ates slightly along the outer wall upstream of the region of turning, but 
unlike example I, a larger deceleration also occurs along the inner wall 
in the downstream half of the elbow. 
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Hie pressure coefficient (p - P-j_) along the irmar »wd outer walls 
far examp le II and for the potential flow solution are presented In fig- 
ure 13 • The pressure coefficient for example H converges to zero far- 
ther upstream of the region of turning than for potential flow. For 
reasons already discussed under example I the exit" pressure' coefficient 
is greater for example II than for the potential flow solution. 


SUPERPOSITION OF SOHJTIONS 

The solutions for example I and for potential flow can he combined 
in linear combinations using the principle of superposition of solutions. 
In example II the principle .of superposition of solutions does not apply, 
because for this solution the right-hand side of the flow equation (10) 
is a function of the dependent variable V* 

Pure shear flow for example I . - The streamlines for pure s hear 
flow are obtained by subtracting the stream, function for potential flow 
from that for example I. Thus the pure shear flow solution is the solu- 
tion for which the net mass flow from wall to wall across lines of con- 
stant £ ls equal to zero. The stream function for pure shear flow is 
zero on both walls and negative throughout the channel. Lines of con- 
stant percentage of minimum stream function are presented in figure 14. 
The inlet velocity profile for the pure shear flow is linear varying 
from -0.25 to 0.25. Velocities are negative along the inner wall and 
positive along the outer wall. The average velocity across any Una of 
constant £ is zero and hence the flow is an eddy occupying the entire 
channe l. The streamlines all begin and end. at the inlet; the zero 
strea m li ne s (walls) extend infinitely far downstream. The pure shear 
flow has no Important physical significance and is of interest only 
mathematically. The solution corresponds to fluid motion in the elbow 
with zero average flow rate but with the elbow rotating in its plane 
with a constant angular velocity equal to one: half the vartlcity of 
example I . __ 

S eparated flow . - If 25 percent, for example, of the potential flow’ 
solution is added to the pure shear flow just discussed, the case of a 
separated flow on the inner wall upstream of the region of turning ls 
obtained. The streamlines for this flow are presented in figure 15. The 
separated flow upstream of. the region of tur ning is forced to reattach by 
the acceleration of the flow -through the elbow. The solution can also be 
reversed, in which case the flow enters at the end with smaller cross 
section; and the solution in figure 15 shows flow separation resulting 
from deceleration. Thus, some aspects of boundary-layer behavior can be 
d em o ns trated by shear flow' solutions that ignore viscous forces in the 
immediate vicinity of the solution. 
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Various percentages of pure shear flow . - The potential flow solu- 
tion was added to various percentages of the pure shear flow solution, 
resulting in the Inlet velocity profiles shown In figure 16. These 
velocity profiles are for flows with varying amounts of vorticity. The 
velocities along the inner and outer walls are shown in figures 17(a) 
and 17(h) . The velocities decrease on the inner wall a nd increase an 
the outer wall as the negative vorticity increases so that the curves 
far the inner and outer walls cross, negative values of Q along the 
i nne r wall would indicate the presence of a separated flow region on 
this wall. The deceleration before turning becomes more pronounced as 
the negative vorticity increases. 

The pressure coefficient (P - p.j_) along the channel walls 1 b 
shown in figure IB. The pressure coefficient along each wall Increases 
as the vorticity increases. As a result, the static pressure drop 
through the elbow becomes less as the percentage of pure rotational flow 
becomes larger. The reasons far this behavior have already been dis- 
cussed under example I. 

The difference in pressure coefficient across the p-hann^i at any 
line of constant £ is nearly the same for all the cases considered in 
figure IS. This fact is shown in figure 19 where the difference in 
pressure coefficient A(P - P^) is plotted against £ . 


SUMMARY QF EESUITB 

A method of a n alysis is developed far two-dimensional shear flows 
in channels of arbitrary shape. Solutions were obtained far shear flows 
in a 90° elbow with linear and sinusoidal velocity distributions across 
the inlet. Several 1 1 near combinations of the potential flow solution 
and the solution for 1 inear inlet velocity distribution were maflp to 
obtain: (l) the pure shear flow solution (negative vorticity), (2) a 

solution that is representative in same respects of the behavior of a 
separated boundary layer under accelerating flow conditions, and (3) a 
series of solutions with equal Increments of negative vorticity. Tht* 
solutions indicated that: 

1. The drop in static pressure across the accelerating elbow 
decreased as the vorticity of the shear flow increased. This behavior 
is explained by the increased nonunif ormity of the inl et velocity with 
increased vorticity. 

2. For linear inlet velocity distributions local deceleration along 
the outer channel wall increased as the negative vorticity of the s hear 
flow increased. However, for the sinusoidal velocity distribution local 
decelerations occurred on both walls. 
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3. Par linear inlet velocity distributions , as the negative vartic- 
ity of the shear flow Increased, the magnitude of the velocities 
increased on the outer wall and decreased on“the inne-r wall, the local 
pressure coefficients increased gradually, and the difference in local 
pressure coefficient on the two walls at corresponding points was 
practically unchanged. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 16, 3S52 
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APPENDIX A 
8YMB0U3 

The follow ing symbols axe used In this report: 

A * area of region in plane of elbow 

g acceleration due to gravity 

H Bernoulli constant, dimensionless, equation (3a) 

h Bernoulli constant, dimensional, equation (3) 

N distance along outer normal to streamline, dimensionless, 

expressed as ratio of v e 

P static pressure, dimensionless, equation (lc) 

P-P^ local pressure coefficient, equation (18) 

AtP-Pi) difference in local pressure coefficient across channel along 
line of constant | 

p static pressure, dimensional 

Q resultant velocity, dimensionless, equation (3b) 

q resultant velocity, dimensional 

S distance along streamline, dimensionless, expressed as ratio 

of w e 

U velocity component in X-direction, dimensionless, equation (Id) 

u velocity component in X-direction, dimensional 

V velocity component in Y-direction, dimensionless, equation (id) 

v velocity component in Y-direction, dimensional 

v width of elbow, dimensional 

X,Y Cartesian nnrrr rH na tes, dimens io nle ss, equation (id) 

x,y Cartesian coordinates, dimensional 
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i vorticity, equation (11) 

r\ stream function for potential flow solution, equation (12a) 

9 angle "between X-axis and streamline measured positive In 

counterclockwise direction 

t velocity potential function for potential flow solution, 

equation (12a) 

p fluid weight density 

f stream function for rotational flawy dimensionless, 

equation (6a) 

Subscripts: 

0 conditions along tj ■ 0 

av average condition across the channel 

e conditions at exit, downstream at infinity 

1 conditions at inlet , upstream at infinity 

mi n minimum 

o conditions in potential flow solution 

component in tj -direction — - - 

conqponent in | -direction 


T} 

e 
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APPENDIX B 


TRANSFORMATION OF DIFFERENTIAL EQUATION OF FLOW FROM 

XT- TO -PLANE 

By def init ion. I and tj satisfy the Cauchy-Riemann differential 
equations, equations (12a) , so that 


£i + sii-o 
ax 2 3 y 2 


afn + afa 

3x z ay 2 


o 


> 


(Bl) 


In order to transform the coordinates of equation (10), the following 
relations were used 


a 2 ? 

ax 2 


d|2 


+ 2 


a 2 ¥ a* a* 

SfSn SxSx 


*3(s)‘ 


arafi arafa 
si ax 2 ax 2 


= IL 


2 a£f 
st 2 


2U 0 V 0 


a 2 ¥ . v 2 a 2 !? . a^r a 2 ! . a? a^i 
° an 2 ^ ax 2 ^ ax 2 


(B2a) 


a% _ a^y /a A 2 
ax 2 ” ai 2 Vs?/ 


+ 2 


a 2 ? ai an ^ a 2 ur /an \ 2 . a 2 e a? . a 2 n anr 

Sf5n 5 y Sy + ^2 \Sy/ + St 2 SI + aT 2 Sn 


= Vo 2 + ZtJoYo 


|!^ + u 0 2 i% + ||^i + p^T L 

a^ 2 ay 2 °ti ay 2 


(B2b) 


From equations (10), (Bl), (B2a), and (B2h) 
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where 



Equation (13) Is the differential equation of flow in |,T) -coordinates . 
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APPENDIX C 


DERIVATION OF EQUATIONS FOR RESULTANT VELOCITY Q, AND ITS 
COMPONENTS Qg, U, AND V 

The derivatives of IF with respect to £ and tj expressed in 
terms of X and Y are given by 


d¥ cfg ax ^ 
Sl = Sx Sf + 


and 


S¥ _ S¥ 5X . 
Stf Sx St[ + 


But, it can be shown that 


3X °°° 

51' So 

5j sin e 0 
5f — 

ax - sln e o 
SY _ cos 3 q 

^ = Qq 


so that from equations (6a), (Cl), and (C2) 


>vr cos 0 o sin 0 O 

5».. T _o +D _o 


d¥ „ 8111 e o . „ cos e o | 

5? ' v ST~ J 


d¥ SY 
Sy SI 




S¥ Sy 
Sy Stf 






(Cl) 


(C2) 


(C3) 
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Pram figure 1 it is seen that 
ft | - Q cos(0-0 Q ) 

■ ft(cos 0 cos 0 Q + sin 6T sin 0 O ) 
= U cos 0 Q + V sin 0 Q 
= ft sin (0-0 Q ) 

=> ft(sin 0 cos 6 q - sin 0 O cos 0) 
= V cos 0 O - U sin 0 O 


Finally, from equations (C3) and (C4) 




(14b) 


^ ■ ■*» If 

(14c) 

Likewise, from equation (C4) 


* 

U = ftj 
V = ft| 

X 

cos 0 O - Q_ sin 0 O 

. 

sin 0 O + ft^ cos 0 O J 

(15) 
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APPENDIX D 


INLET AND EXIT CONDITIONS FOR EXAMPLE I 


At the inlet of the elbow, upstream, at infinity, Qq is l/2, 
V 10 o, H is eqiyal to and 7 is a function of tj only, bo 
that from equations (14b) and (16a) 




(0.25 + 0.50 r^) cLTli 


7 ± * 0.5(t} ± + tj ± 2) 


(I6d) 


Also, at the inlet P is constant and H and Q are functions ' 
of T] only, so that from equation (4) 



and, therefore, from equation (16 a) 



+ 0 . 50 ) 0.5 dT]^ 


% = + 0.25(th + TJi 2 ) 


or from equation (l6d) and the fact that H is a function of “ST only 

H = Hq + I (16b) 

u 

At the. exit of the elbow, downstream at infinity, Qq is 1.0, 

Qri is 0, Qj is equal to Qg, and 7 is a function of t] only. 
Therefore, from equation (13) 
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which, after integration, becomes 

v e 0.875 Tj e - + 0.125 T} e 2 (ige) 

Therefore, from equations (14b) and (l6e) 

Qg = 0.875 + 0.25 rj e ' (l6f) 
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APPENDIX E 


INLET AND EXIT CONDITIONS FOR EXAMPLE II 

/ 

At the inlet of the elbow, upstream, at infinity, Qq is l/2, 

Qij is 0, Q| is equal to Qi, and 7 is a function of tj only, bo 
that from equations (14b) and (19a) 


r i pm 

Jo 


sin nrq ^ dT^ 


V± = h (l-cos rt7] ± ) 


(I9d) 


Also, at the inlet P is constant and H and Q are functions 
of tj only, so that from equation (4) 


20 , (SS\ = [&) 

T- \dri/i 


and, therefore, from equation (19a) 


OH, 


®i -T 


Hi 


Hi,0 


sin JCTJi COS JCTli dT]i 


2 

Hi = Hi^ 0 + sin 2 jnji 


or from equation (l9d) and the fact that H is a function of 7 only 

H = Eq + ^ (T-IF 2 ) (19b) 

At the exit of the elbow, downstream at infinity, Qg is 1.0, 

Oq is 0, Qt is equal to Qg, and 7 is a function of tj only. 
Therefore, from equation (13) 


(: 


dh \ _■£_ jc 2 ¥ 

aWe ' 8 4 
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which, after integration between -limits, becomes 


Tr e = | (i - cos + s in ^n e ) 
Therefore, from equations (14b) and (l9e) 


(I9e) 


Qe - r ( sln f^e + cos fife) 


(I9f) 
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(a) lunar wall. 
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Figure 17 • - Velocities along inner and outer walls for potential flow plus various percentages of pure shear flaw. 





e 

(b) Outwr wall. 



Jignr* 17. - Concluded. Velocities alone timer and outer valla far potential flow plus various percentages of pure shear flow. 
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